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ACC.0002US 

USING MAGNETIC RESONANCE IMAGING TO DIRECTLY MAP NEURONAL 

ACTIVITY 

This application claims priority to United States 
Provisional Application Number 60/412,171 filed on 
September 20, 2 0 03 in the names of Jinhu Xiong, Jia-Hong 
Gao, and Peter T. Fox, entitled "Using Magnetic Resonance 
5 Imaging to Directly Map Neuronal Activity". 

Background 

The present invention relates to medical imaging and 
more specifically magnetic resonance imaging (MRI) . 
Various neuroimaging techniques are presently 

10 available. These techniques include MRI, functional MRI 
(f MRI) , and positron emission tomography (PET) , for 
example. None of these techniques, however, are able to 
directly measure neural activity, i.e., brain activity. 

Instead, these techniques detect brain activity via 

15 cerebral hemodynamic and metabolic responses to neural 

firing. However the temporal resolutions of f MRI and PET 
are ultimately limited by the slow response function of 
cerebral hemodynamics, which is on the order of seconds. 
Furthermore, their inferences regarding neuronal activity 

2 0 are necessarily complicated by the variability of coupling 
between neuronal activity, cerebral hemodynamics, and 
metabolism. 
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Other techniques to map neural activity include 
electroencephalography (EEG) and/or magnetoencephalography 
(MEG) . However, these techniques often have poor spatial 
resolution. Because both EEG and MEG rely on information 
5 detected at the scalp to localize active sites inside the 
brain, both EEG and MEG require solving an inverse problem, 
which leads to spatial uncertainty in the localization of 
electromagnetic sources. In addition, EEG and MEG are each 
limited in the activation geometries they can detect and 

10 are unable to detect neuronal activities deep in the brain. 
While combining information from modalities detecting 
different physiological variables (for example, data from 
f MRI and MEG) can partially offset the drawbacks of the 
individual modalities and can provide brain activation maps 

15 with high spatial and temporal resolution, the basic 

limitations for each modality, such as the indirect nature 
of f MRI measurement and the inverse problems for EEG and 
MEG, remain obstacles. 

Thus a need exists to provide improved combined 

20 spatio-temporal resolution of neural activity imaging. 

Further, a need exists to directly map such neural activity 
to avoid the inverse problem, and directly measure magnetic 
sources originating from neural firing with high spatio- 
temporal resolution . 
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Summary of the Invention 
Functional operations of the human brain are 4- 
dimensional (4D) processes (i.e., time and space). Even 
simple tasks activate neural populations distributed across 
5 the brain in space and time, which can be fully mapped only 
by techniques with high spatio-temporal resolution. 

Embodiments of the present invention include methods 
and apparatus to directly map neural activity using MRI . 
As used herein the terms "directly mapping" or "direct 

10 mapping" mean the measurement of neural activity 

concurrently with neural electromagnetic changes. That is, 
in contrast to other techniques in which such neural 
activity may be inferred or measured via cerebral 
hemodynamic and metabolic changes (incurring a temporal 

15 delay) , "directly mapping" means the measurement of 
electromagnetic effects of neural activity without a 
temporal delay and without a need for solving an inverse 
problem. This direct mapping may be referred to herein as 
"magnetic source MRI" (msMRI) . In so doing, system- level , 

20 event-related neuronal activity may be mapped with high 
spatial and temporal accuracy by directly detecting 
magnetic transients induced by neural firing. Further, via 
such temporal accuracy, latency for activation to spread 
from location to location may be measured. Such latency 

25 may include absolute latency (e.g., relative to a delivered 
stimulus); and inter- regional latency (i.e., signal travel 
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time between brain regions) . Such measurements of 
intracerebral latencies may be used to model neural system 
and may aid in diagnosis and study of brain disorders. 

Embodiments of the present invention may be used to 
5 greatly increase the temporal resolution with which MR I can 
map brain functional activity and provide a powerful new 
tool for mapping brain functional organization in humans 
and animals. 

Further still, embodiments of the present invention 

10 may be used to detect intrinsic rhythms within the nervous 
system. That is, in accordance with the present invention, 
embodiments may detect intrinsic brain oscillations which 
may be present in the absence of perturbations. In so 
doing, changes in intrinsic rhythms over time may be 

15 measured and analyzed. These changes may be used to 

indicate various states of nervous system functionality. 
For example, such changes in rhythm may relate to brain 
damage, drug effects, and or nervous system 
disease/disorders . 

2 0 In other embodiments, similar analyses of nervous 

system state may be accomplished by use of msMRI in 
connection with a mental event. As used herein the term, 
"mental event" means any externally or internally caused 
activation, perception, cognition, emotion and or 

2 5 attention. Further, as used herein the term "event 

related" may mean both block design activities in which 
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multiple stimuli are presented and measured, and event 
related activities in which a stimulus is provided and a 
series of measurements are then made. 

Further still, in certain embodiments, msMRI may be 
5 conjoined with one or more other measurements of nervous 
system state. For example, msMRI may be performed in 
connection with fMRI to measure cerebral hemodynamic, 
metabolism and/or neuronal activity. Similarly, msMRI may 
be performed in connection with EEG or MEG. In yet other 

10 embodiments, msMRI may be combined with other imaging 
techniques, such as PET, computed tomography (CT) , or 
single-photon computed tomography (SPECT) , for example. In 
so doing, metabolic, hemodynamic, and electromagnetic 
activity of the nervous system may be measured. 

15 By using embodiments of the present invention, co- 

variants between different neural sites may be measured and 
analyzed. Such co-variants may be used in diagnosis of 
disease or disorders of the nervous system, as well as the 
analysis of drug effects and analysis of anatomical 

20 connectivity between neural sites. 

In one embodiment, the present invention may include 
detecting, using MR I , regional neural activity in a subject 
undergoing MR I scanning, based on magnetic fields induced 
by the regional neural activity; and spatially and 

25 temporally localizing the neural activity using at least a 
portion of the detected magnetic fields. 



5 



Brief Description of the Drawings 
FIG. 1 is a flow chart of an example method according 
to the present invention. 

FIG. 2 is a block diagram of a MR I system in 
5 accordance with an embodiment of the present invention. 

FIGS. 3A and 3B are graphical representations of 
hemodynamically neutral data acquisitions according to the 
present invent ion . 

FIG. 4 is a time sequence for an exemplary study in 
10 accordance with the present invention. 

FIG. 5 shows spatio-temporal plots of event-related 
neural activity of a subject detected by msMRI and BOLD 
fMRI maps for the subject in accordance with the present 
invention. 

15 FIG. 6 is a graphical representation of averaged time 

courses of relative activation sizes in the exemplary study 

according to the present invention. 

FIG. 7 are functional activation maps and graphical 

representations of event related signal changes for two 
2 0 subjects in the exemplary study according to the present 

invention. 

FIG. 8 is a graphical representation of signal change 
versus echo time in accordance with one embodiment of the 
present invent ion . 
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FIG. 9 is a representation of a two dipole system 
showing parallel (left) and anti-parallel (right) 
configurations . 

FIG. 10 is a representation of two adjacent two dipole 
5 systems in an anti-parallel configuration. 

FIG. 11 is a graphical representation of a percentage 
of signal change versus number of dendrites for parallel 
and anti-parallel configurations in accordance with one 
embodiment of the present invention. 
10 FIG. 12 is a graphical representation of computed MR I 

signal changes with echo time for parallel and anti- 
parallel configurations in accordance with an embodiment of 
the present invention. 

Detailed Description 

15 In one embodiment, mapping brain activity with MR I by 

detecting magnetic fields induced by neural firing may be 
accomplished by inducing and detecting phase coherent 
signals of proton nuclear spins. Neuronal activity creates 
ionic currents which induce weak electromagnetic fields 

20 (~10~ 13 Tesla (T) ) . Nuclear spins exposed to neuronally 
induced magnetic fields will lose phase coherence, which 
will slightly decrease MR I signal strength. Neuronal 
magnetic transients may be mapped, therefore, by detecting 
event-related decrements in the MR I signal. 

2 5 Referring now to FIG. 1, shown is a flow chart of an 

example method according to the present invention. As 
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shown in FIG. 1, the method begins with the location or 
placement of a subject within an MRI scanner (block 100) . 
It is to be understood that any MRI scanner may be used in 
connection with the present invention. As is well known in 
5 the art, such subject location may place the subject in an 
RF coil. Furthermore, the MRI magnets may be shimmed to 
the subject. Furthermore, within block 100 it is to be 
understood that other well known techniques in the 
initialization of MRI scanning may be performed, such as 

10 preparing a scout image. 

Next, the subject is stimulated and data is obtained 
(block 110) . While not shown in FIG. 1, it is to be 
understood that obtaining data and stimulating the subject 
may be performed iteratively, as desired by a particular 

15 msMRI procedure. While it is to be understood that in 
certain embodiments, the stimulation may be an external 
stimuli such as a task (e.g., visuomotor task), in other 
embodiments, the stimuli may be internally provided. 
Further still, in other embodiments, no stimuli need be 

20 presented, and msMRI may measure intrinsic rhythms of the 

nervous system in the absence of stimuli. Data is obtained 
in accordance with well known MRI techniques. The data 
obtained may include, for example, information regarding 
transient magnetic fields. Information of transient 

25 magnetic fields may be obtained by measuring MRI signal 

changes in magnitude or phase. In certain embodiments the 



data obtained may relate to MR I signal magnitude. Because 
orientations of neuronal currents are very complex and 
pseudo- random, phases may destructively add, resulting in a 
minimal phase shift and a large phase dispersion. Magnitude 
5 measurements, therefore, may be more sensitive for 

measuring the collective neuronal activity than phase 
measurements . 

Finally, data processing is performed to obtain the 
desired msMRI images (block 120) . Such processing may be 

10 performed in accordance with well known MR I techniques in 

which images are reconstructed from the data obtained using 
a data analysis machine. In certain embodiments the data 
analysis machine may be data processing system, such as a 
personal computer, workstation, or the like. The data 

15 processing typically includes multiple steps, for example, 
image reconstruction, motion correction, spatial 
normalization, value normalization/correction, statistical 
analysis, and statistical inference. 

Various data processing techniques may be used to 

2 0 analyze the images. For example, various image 

reconstruction methods may be performed to obtain desired 
images from the MRI signal. Various motion corrections may 
be performed to minimize image artifacts related to subject 
movements. Various spatial normalization procedures may be 

25 performed to standardize the size, position, and 

orientation of images across different subjects. Various 
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value normalization/correction methods may be performed to 
minimize the physical and physiological artifacts of the 
MRI signal acquired over time. Various statistical 
analyses may be performed to create different statistical 
parametric images such as z- , t-, or r- images. 

In certain embodiments, the data processing techniques 
may be implemented in software. As such, these embodiments 
may be stored on a storage medium having stored thereon 
instructions which can be used to program a data processing 
device, such as a computer system or the like, to perform 
the embodiments. The storage medium may include, but is not 
limited to, any type of disk including floppy disks, 
optical disks, CD-ROMs, CD-RWs, and magneto-optical disks, 
semiconductor devices such as ROMs, RAMs, EPROMs, EEPROMs, 
magnetic or optical cards, or any type of media suitable 
for storing electronic instructions. Similarly, embodiments 
may be implemented as software modules executed by a 
programmable control device. A programmable control device 
may be a computer processor or a custom designed state 
machine, for example. 

It is to be understood that various MRI pulse 
sequences may be used. In various embodiments of the 
present invention, an asymmetric pulse sequence may be 
used. Such an asymmetric pulse sequence may be used to 
avoid cancellation of the signal by the second pulse in the 
pair. In one embodiment, a gradient -echo echo-planar- image 



10 



(EPI) pulse sequence may be used. In various embodiments, 
the asymmetric pulse sequence may have a repetition time of 
between approximately 40 and 10,000 milliseconds, an echo 
time of between approximately 10 and 200 milliseconds, and 
5 a flip angle of between approximately 10 and 180 degrees. 
More specifically, in certain embodiments, the pulse 
sequence may include a repetition time of between 
approximately 4 00 and 4,000 milliseconds, an echo time of 
between approximately 50 to 150 milliseconds, and a flip 

10 angle between approximately 40 to 120 degrees. However in 
other embodiments, other parameters may be used. Further 
in other embodiments, other pulse sequences such as a 
conventional spin echo, gradient echo, fast spin echo, 
spin-echo EPI, or a spiral imaging pulse sequence may be 

15 used. 

In certain embodiments, a hemodynamically neutral 
msMRI technique may be used. Brain activity induces large 
hemodynamic effects (e.g., increased blood flow). These 
hemodynamic effects are slow but strong. Generally, they 

20 reach a peak only after 4 to 6 seconds and decline equally 
as slowly. Thus as used herein the term "hemodynamically 
neutral" means the obtaining of data during a time in which 
hemodynamic effects are at a substantial steady state. 

Referring now to FIG. 2, shown is a block diagram of 

25 an MR I system in accordance with an embodiment of the 

present invention. As shown in FIG. 2, system 200 includes 
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an MRI scanner 2 05 having a gradient coil 210 and an RF 
coil 220. Shown in FIG. 2, the gradient coil 210 may be 
controlled via x, y, z amplifier 215, while the RF coil 220 
may be controlled via an RF amplifier 225. While shown as 
5 including such coils, it is to be understood that FIG. 2 is 
shown for illustrative purposes and an MRI scanner in 
accordance with an embodiment of the present invention may 
include additional components of a standard scanner such as 
a main magnet, additional gradient magnets, and the like. 

10 The amplified signals provided by MRI scanner 2 05 may 

be provided under control of a controller 250. In one 
embodiment, controller 250 may be a workstation, server, 
personal computer or other data processing system capable 
of providing control signals for system 200. 

15 Coupled to controller 250 may be a display 260 and an 

input device 2 65 for providing display information 
regarding an MRI process and to receive input information 
from a user, respectively. More so, a receiver 240 may be 
coupled to MRI scanner 2 05 to receive resonance signals 

2 0 therefrom and preprocess them. As shown in FIG. 2, 

received resonance signals may be provided from receiver 
240 to controller 250 for further processing. Alternately, 
such resonance signals may be provided to a separate data 
processor for desired processing. 

2 5 In accordance with an embodiment of the present 

invention, controller 250 may include software routines and 



the like to obtain msMRI signals and use them to map neural 
activity with a high degree of spatial and temporal 
localization. 

While shown with the specific components discussed 
5 above in FIG. 2, it is to be understood that embodiments of 
the present invention may be used with various MR I devices 
now known or available in the future. 

Referring now to FIGS. 3A and 3B, shown are graphical 
representations of two manners of obtaining a 

10 hemodynamically neutral response in accordance with the 
present invention. As shown in FIG. 3A, a blood oxygen 
level -dependent (BOLD) signal 10 rises to a peak 
approximately five seconds after an event. To provide for 
a hemodynamically neutral analysis, data may be obtained in 

15 windows 15 and 20, occurring prior to and after hemodynamic 
activity. Alternately, in an msMRI technique in which 
multiple rapid stimuli are provided, the hemodynamic 
response causes fluxes in brain activity too fast for 
hemodynamic changes to track. That is, a hemodynamic 

2 0 response plateaus and remains at a steady state throughout 
the data acquisition. In such a technique, a BOLD signal 
3 0 remains at a substantial steady state, so that data 
acquisition may occur in window 35, shown in FIG. 3B. 

Exemplary Study 

25 In a study according to the present invention using a 

well-established visuomotor paradigm, human brain msMRI 
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images correctly detected the locations and latencies of 
activations in visual, motor, and premotor cortices, with a 
temporal resolution of 100 milliseconds (ms) and spatial 
resolution of 3 millimeters (mm) . Signal strength was 
5 comparable to other event-related functional MR I methods: 
about 1% of the baseline signal. 

In this exemplary study, a simple, well-established 
visuomotor task was used to map the system- level 
organization of the human motor and visual cortices using 

10 msMRI . In the study, the subject's heads were immobilized 
in a closely fitted, thermally molded, plastic facial mask 
individually formed for each subject. Such masks minimize 
head movement during MR I scanning. Cued by a brief (50 ms) 
stimulus (a wedge of random dots) in the lower left visual 

15 field, subjects pressed and released a button with the 
right index finger. Magnetic source MR I images were 
acquired over a 1300 ms period starting at 200 ms before 
cue onset and divided into thirteen consecutive, 100 ms 
time frames. Six time frames (0 - 600 ms after cue onset) 

2 0 were reported. Images were spatially configured as five 

contiguous axial slices with a slice thickness of 6 mm and 
orientated to cover both the upper bank of the calcarine 
cortex (the lower field representation of primary visual 
cortex) and the hand area of primary motor cortex. In plane 

25 spatial resolution was 3 mm. A relatively long echo time 
(TE = 100 ms) was used to maximize the strength of the 



msMRI signal. Six normal subjects (five men and one 
woman) were scanned. Three of them were scanned twice one 
week apart to assess reproducibility of brain activation. 
In addition to msMRI, a high resolution (lxlxl mm) 
5 anatomical MR I and a "traditional' 7 BOLD fMRI were acquired 
in each subject. The BOLD fMRI was acquired in a block 
design, using the same visuomotor task and a resting-state 
control . 

Referring now to FIG. 4, shown is a time sequence used 

10 for the exemplary study. For each stimulation ON/OFF cycle 
50, one five-slice image was acquired corresponding to the 
ON period (solid arrows) and another corresponding to the 
OFF period (dash arrows) . Each slice was color-coded. The 
time between adjacent slices was 200 ms . The MRI scanner 

15 was precisely synchronized with stimulation onset. 

As shown in FIG. 4, data acquisitions 60 for the first 
five runs start at -200, -100, 0, 100, and 200 ms relative 
to stimulation onset. In this exemplary study, MRI data 
were acquired on a 1 . 9 T GE/Elscint PRESTIGE whole-body MRI 

20 scanner using a gradient-echo echo-planar- image (EPI) pulse 
sequence with the following parameters: repetition time, 
TR, of 1000 ms, echo time, TE, of 100 ms, and a flip angle 
of 90 degrees. 

Five contiguous, oblique slices were acquired, with an 

25 in-plane spatial resolution of 3 X 3 mm 2 and slice thickness 
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of 6 mm. The orientation and location of the slices were 
carefully selected to include both visual and motor areas. 

For the exemplary study, each MR I session included six 
runs of data acquisitions, with an acquisition time of 200 
5 seconds per run. The first five runs were designed for 
msMRI imaging of event -related neuronal activity and are 
shown in FIG. 4. Each run consisted of 100 ON/OFF cycles, 
with two multiple-slice images for each cycle (one inter- 
stimulation interval (ISI) ) . The last run was a control. 

10 It was acquired either at a resting-state (subject does not 
actively perform any task) or as a block design 
conventional blood oxygen level dependent (BOLD) f MRI 
study, in which subjects performed the visuomotor task for 
the first 90 seconds and then rested for the rest of time. 

15 Referring now to FIG. 5 (left side) , shown are spatio- 

temporal plots of event-related neuronal activity detected 
by msMRI for a single subject in accordance with the 
exemplary study. For comparison, the same subject's BOLD 
f MRI maps are also shown in FIG. 5 (right side) . The time 

2 0 sequence of visual and motor events is shown along the top 
of FIG. 5. The activation information is overlaid on Tl- 
weighted MRI images acquired at the same location and 
orientation. Because different MRI slices were acquired at 
slightly different times, spatial normalization was not 

25 performed to preserve temporal information. The letter L 
on the left-lower corner indicates the left cerebral 
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hemisphere. The letter P indicates posterior. The term VI 
refers to the right striate cortex, SI refers to the left 
somatosensory cortex, SMA refers to the anterior 
supplementary motor area, and Ml refers to the left primary 
5 motor and premotor cortices. While not shown in FIG. 5, 
color images may be used in which a color scale exists to 
represent the t value of each voxel . 

For the exemplary study, the MR I images were processed 
using image processing software. Several modifications to 
10 the software were performed to accomplish this study. For 
example, motion correction, and value 

normalization/correction procedures were modified, as 
discussed above. 

In the exemplary study, the first 2 0 images of each 

15 run were discarded to allow hemodynamics and the MRI signal 
to reach a steady state. All data were assessed for inter- 
scan, intra-subject movement. A two-dimensional (2D) 
movement correction was performed to minimize in-plane 
motion. Data interpolation between image slices was 

20 purposely avoided because different slices were acquired at 
different times. Voxel -by -voxel linear detrending was 
performed to remove the linear drift of the MRI signal. A 
mean image was created for each on-off cycle by averaging 
across the time-series in the cycle; the mean image was 

25 then subtracted from each image to create residual images. 
A 2D spatial Gaussian filter with a full width at half 
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magnitude of 4.5 mm was applied. A group Student's t test 
was performed on the residual images. The t-image was then 
thresholded using a t -value threshold of 3 . 0 (p < 0.0013) 
and cluster size threshold of 4 voxels to detect 
significant activation . 

Referring now to FIG. 6, shown is a graphical 
representation of averaged time-courses of relative 
activation sizes for different brain regions in accordance 
with the exemplary study. Activation sizes were averaged 
across subjects and sessions and were normalized by 
dividing each size with the maximum size for the region. 
Four different brain regions are shown: the right striate 
cortex (VI) , the left primary motor and premotor cortices 
(Ml) , the left somatosensory cortex (SI) , and the anterior 
supplementary motor area (SMA) . Time frames 1 -5 represent 
time intervals of 0 - 100 ms, 100 - 200 ms, 200 - 300 ms, 
300 - 400 ms, and 400 - 500 ms after the visual stimulation 
onset , respectively . 

Referring now to FIG. 7, shown are functional 
activation maps and graphical representations of event - 
related signal changes for two different subjects in 
accordance with the exemplary study. As shown in FIG. 7, 
the activation information has been overlaid on Tl-weighted 
MR I images acquired at the same location and orientation. 
Again, while not shown in FIG. 7, for color images, a color 
scale may be used to represent the t value of each voxel . 
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The images shown in FIG. 7 are oriented the same as those 
of FIG. 5. The graphical representations represent signal 
changes corresponding to the activated areas inside the 
boxes. Each data point represents an average of 18 
individual trials . 

A current -dipole model was constructed for modeling MR I 
signal changes resulting from neuromagnetic fields. The 
model was constructed over a scale range from a single 
neuron to a typical MR I voxel (approximately one million 
neurons) . Each dendrite was modeled as a current dipole. 
Distributions of neuronal magnetic fields were estimated and 
the interaction between the neuronal magnetic fields and 
nuclear spins was assessed across a range of scales, 
orientations, configurations, and distributions of dendrite 
packing density. More detailed information regarding 
modeling is discussed below. 

As predicted, msMRI detected regional signal 
decrements in visual, sensorimotor, and pre-motor cortices 
with appropriate latencies (below), locations, and 
lateralities. For all subjects, occipital activations 
bordered the calcarine fissure and were predominately right 
hemispheric, consistent with the left visual field location 
of the cue. Similarly, sensorimotor and premotor 
activations were chiefly left hemispheric, consistent with 
the right-hand motor response. Primary motor area (Ml) 
responses lay immediately anterior to the central sulcus; 
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primary sensory cortex (SI) responses lay immediately 
posterior to the central sulcus. The locations of task- 
induced activations detected by msMRI and BOLD fMRI were in 
good agreement in all subjects. In the illustrated subject 
5 of FIG. 5, locations and lateralities of msMRI activations 
in Ml, SI, the supplementary motor area (SMA) and posterior 
cingulate were closely replicated by BOLD fMRI. 

In all areas imaged, msMRI signals were detected 
predominately in cerebral grey matter, as shown in FIG. 5. 

10 Grey matter localization of msMRI signals may appear 

counter- intuitive, but is in good agreement with known 
electrophysiology . 

In addition to accurate spatial information, msMRI 
maps provided richly detailed temporal information 

15 regarding task- induced neuronal activity. As illustrated 
in FIG. 5 (single subject) and FIG. 6 (group data) , msMRI 
maps showed right-hemispheric activation of primary visual 
(striate) cortex (VI) in frame 1 (0 - 100 ms) immediately 
after visual stimulation onset (FIG. 5) . The visual 

20 activation moved laterally to the right extrastriate visual 
areas (Brodmann area [BA] 18 and 19) in frame 2 (100-200 
ms) . Both the striate and extrastriate visual cortical 
areas were re-activated in frame 4 (300 - 400 ms) . Onset 
latencies of msMRI visual activations agree well with known 

25 EEG onset latencies. Re-entrant activation of visual cortex 
has also been reported, with onset at 250 - 300 ms, about 
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50 ms earlier than the present study. Note that msMRI 
measures integrated activation within a frame (100 ms) . 
Activation starting late in a frame may not be detected 
until the subsequent frame, which may explain why 
5 activation onsets from EEG are sometimes shorter than those 
from msMRI . The visual activation patterns illustrated in 
FIG. 5 were fairly consistent across subjects and sessions. 
Group data clearly showed the same activation-reactivation 
patterns with the similar latencies (FIG. 6) . 

10 In the motor system as well, msMRI detected activation 

spatially and temporally discrete and appropriate to the 
task. In the illustrated subject of FIG. 5, the left 
primary motor (Ml) and premotor areas (BA 6) were activated 
twice: first in frame 2 (100 - 200 ms) and then in frame 4 

15 (300 - 400 ms) . Similar activation patterns and onset 

latencies were observed in the group data (FIG. 6) . The 
first Ml activation is likely preparation for and execution 
of the button press. Using a similar visuomotor task, onset 
latencies of 130 - 180 ms for motor and premotor cortices 

2 0 have been reported, in excellent agreement with the present 
study. The second Ml activation likely represents the same 
components for button release. Chronometric measurements 
showed an average delay of 185 ms between button press and 
button release, which closely matches the delay between the 

25 first and second Ml activations. As expected, the onset of 
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brain activation preceded finger movement by about 100 ms 
(finger response time: 239 ± 57 ms) . 

As shown in FIG. 5, primary somatosensory (SI) cortex 
activation (0 - 100 ms) preceded Ml activation, and 
5 continued during Ml activation (100 - 2 00 ms) , but reached 
maximum after Ml activation (200 - 300 ms) . Similar 
activation patterns and onset latencies were observed in 
the group data (FIG. 6) . Judged by latency, the strongest 
SI activation likely reflects the sensory feedback from the 

10 finger movement. Left SI was activated again during the 
button release period (at 300 - 400 ms) (FIG. 5) . Also 
activated significantly was the anterior supplementary 
motor area (SMA) (BA 6) , which activated first in frame 1 
(0 - 100 ms) , reached maximum in frame 2 (100 - 200 ms) , 

15 and continuously activated in frame 3 (200- 300 ms), frame 
4 (300 - 400 ms) , and frame 5 (400 - 500 ms) (FIGS. 5 and 
6) . The posterior cingulate (BA 31) was activated first in 
frame 1 (0 - 100 ms) and again in frame 3 (300 - 400 ms) 
(FIG. 5) . These cingulate activations likely correspond to 

20 the Bereitschaf tpotential . 

The consistency of msMRI in detecting neuronal 
activity was assessed across subjects and sessions. For all 
subjects, background noise level in the control studies 
(the resting state) was very low and no false-positive 

25 activations were detected. Trial-by-trial consistency was 
excellent (FIG. 7) , with a clear distinction between 
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baseline and activation. Visual inspections of activation 
maps showed good consistency across subjects and sessions. 
For example, for the two subjects shown in FIG. 7, the 
locations of activations in Ml, SI, premotor, SMA, and 
5 posterior cingulate cortices are quite similar. In this 
study, quantitative assessments of activation location 
variability were not performed, as this requires a much 
larger sample size and requires spatial normalization. 

Data according to the present study contradict the 

10 prediction that the magnitude of msMRI signals in the human 
brain would be very weak. Specifically, in the present 
study, msMRI signals are similar in magnitude to event- 
related BOLD fMRI (-1%) . For example, the average 
magnitude (n = 6) of the msMRI signal in left Ml cortex was 

15 1.12% ± 0.54% of the background anatomical signal (FIG. 7). 
At the same field strength (1.9 T) , event-related BOLD fMRI 
effects have been reported at -1% of the background signal. 
The BOLD fMRI protocol, however, used an echo time (TE) of 
4 5 ms, as compared to 10 0 ms TE for the present study. 

2 0 Functional MRI signal changes are TE dependent, with longer 
TE yielding larger percentage signal changes. 

In various embodiments in accordance with the present 
invention, a long TE may provide for increased magnitude of 
msMRI signal strength. In certain embodiments, the TE may 

2 5 be between approximately 10 to 200 ms, and more 
specifically between approximately 50 to 150 ms . 
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Further contributing to increased magnitude msMRI 
signals in certain embodiments may be the measurement of 
MR I signal magnitude change rather than MR I signal phase 
change. Because orientations of dendrites are very complex 
5 and pseudo-random, phases may destructively add, resulting 
in a minimal phase shift and a large phase dispersion. Thus 
in certain embodiments, magnitude measurements may be used. 
In such embodiment, these measurements may be more 
sensitive for measuring the collective neuronal activity 

10 than phase measurements. 

The neuronal activity which induces transient magnetic 
fields detectable by msMRI may also cause cerebral 
hemodynamical changes in the surrounding brain tissues. In 
various embodiments, effects of cerebral hemodynamics may 

15 be sought to be minimized, while maximizing the effects of 
the transient magnetic fields. In the present study, the 
stimulation ON/ OFF cycle was switched rapidly so that the 
cerebral hemodynamical response was effectively in steady 
state. The hemodynamical ly based BOLD signal, therefore, 

2 0 was effectively constant during stimulation ON and OFF 
phases . 

In the present study, the temporal resolution of msMRI 
is equal to the echo time (TE) , 100 ms for the present 
study. This temporal resolution is desired, but not 
25 optimal, to investigate neuronal activity at the system 

level. Activation of a neural population generally lasts 
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tens to hundreds of milliseconds. Similarly, latencies of 
activations are also in a range of tens to hundreds of 
milliseconds. Thus, improving the temporal resolution of 
msMRI may be desired in certain embodiments. The temporal 
5 resolution of msMRI is mainly limited by contrast-to-noise 
ratio, which is affected by MR I parameters including field 
strength, repetition time (TR) , and echo time (TE) , as well 
as by the pattern of neural firing. Temporal resolution 
may be further improved by optimization of experimental 

10 design and MR I pulse sequences, in certain embodiments. 

To confirm that the signals detected by embodiments of 
the present invention are the direct effects of neuronal 
magnetic fields, theoretical modeling was performed and 
experimentally confirmed. A current -dipole model, commonly 

15 used in MEG for the estimation of the magnetic fields 

induced by neuronal firing, was constructed over a scale 
range from a single neuron to a typical MR I voxel 
(approximately one million neurons) . In the range of 
parameters proper for human cortex, modeling showed that 

2 0 changes in the magnitude of the MR I signal due to neuronal 
activity could reach a few percent (0.5% to 5%) . That is, 
msMRI signals in accordance with one embodiment of the 
present invention are within a measurable range, comparable 
to standard f MRI techniques . 

2 5 Referring now to FIG. 8, shown is a graphical 

representation of percentage of signal change versus echo 
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time. The solid line represents theoretical results, which 
will be discussed below. The points on the vertical lines 
represent experimental data, which represent the msMRI 
signals averaged across all significantly activated legions 
5 and across a group of three subjects, where the arrow bars 
represent one standard error. In the study, the parameters 
used for calculating the MR I signal changes are: radius of 
the dendrites a = 0.5 \xm, electrical conductivity of 
axoplasma <j± = 2 Q^m" 1 , and membrane potential Av± = 75 mV. 
10 The number of dendrites fired simultaneously at any time 

was presumed to be approximately 0.2 million for a typical 
MR I voxel. This theoretical prediction was confirmed by 
the experimental data, which showed a signal change of 
1.12% ± 0.54% of the background anatomical signal in the 
15 left Ml cortex. Response magnitudes for other brain 
regions were similar. 

Modeling also tested the relationship between the 
msMRI signal strength and the echo time (TE) . A 
counterintuitive prediction of the model is that a non- 
20 linear relationship exists between TE and msMRI magnitude: 
msMRI signal increased by a factor of 3 . 8 when TE doubled. 
This is in sharp contrast to a nearly linear relationship 
between TE and BOLD f MRI signal. The predicted non- 
linearity was experimentally confirmed, as shown in FIG. 8, 
25 which is a reason why a long TE for msMRI may be desired. 
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For neuronal magnetic fields, the direction of the 
field on the left side of a current dipole is always 
opposite to that on the right. A typical MR I voxel 
contains approximately one million dipoles. Both positive 
5 and negative phases of MR I signals will destructively add, 
lowering net phase for an imaging voxel. Thus, the net 
phase is near zero and is undetectable. On the other hand, 
magnitude of MR I signal significantly changes due to this 
neuromagnetic field. Even though millions of dendrites are 

10 synchronously activated, the combined magnetic field 

remains highly inhomogenous . Nuclear spins experiencing 
these local field inhomogeneities will lose phase 
coherence, resulting in a decrease of MR I signal magnitude. 
Magnitude measurements used in accordance with an 

15 embodiment of the present invention may thus be far more 
sensitive for measuring in vivo neuronal activity than 
phase measurements. In addition to electromagnetic 
effects, neuronal activity also induces physiological 
(e.g., changes in blood flow and metabolic rate) and 

20 mechanical (e.g., cell swelling) effects. 

Thus, spatially and temporally precise signals 
detected using a model -based experimental design optimized 
to detect neuronal magnetic effects are, in fact, msMRI 
effects. The msMRI activations showed a spatial and 

25 temporal distribution appropriate to the neural systems 

activated by the widely used visuomotor task. The sign and 
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magnitude of the observed signals (about -1% of baseline) 
were as predicted by the model. The nonlinear effect of TE 
on the observed signal was as predicted by the model . 
Further, non-magnetic effects of neuronal activity have 
5 been addressed but found wanting as explanations of the 
observed effect. 

Techniques in accordance with an embodiment of the 
present invention offer several advantages over current 
neuroimaging methods. It appears to provide better 

10 combined spatio-temporal resolution than any currently used 
non- invasive neuroimaging methods. Compared with 
"traditional" f MRI and positron-emission tomography (PET) , 
msMRI offers much higher temporal resolution, with no loss 
of spatial resolution. As discussed above, detecting brain 

15 activity via the cerebral hemodynamic and metabolic 

responses to neural firing, the temporal resolutions of 
f MRI and PET are ultimately limited by the slow response 
function of cerebral hemodynamics, which is on the order of 
seconds. Furthermore, their inferences regarding neuronal 

20 activity are necessarily complicated by the variability of 
coupling between neuronal activity, cerebral hemodynamics, 
and metabolism. Compared with EEG and MEG, msMRI offers 
higher spatial accuracy. Relying on information detected at 
the scalp to localize active sites inside the brain, both 

25 EEG and MEG require solving an inverse problem, which leads 



to spatial uncertainty in the localization of 
electromagnetic sources . 

In contrast msMRI effects in accordance with one 
embodiment of the present invention are spatially mapped in 
5 the same manner as traditional MR I techniques and involve 
no inverse problem. In addition, EEG and MEG are each 
limited in the activation geometries they can detect and 
are unable to detect neuronal activities deep in the brain; 
msMRI has no such limitation. Combining information from 

10 modalities detecting different physiological variables (for 
example, data from f MRI and MEG) can partially offset the 
drawbacks of the individual modalities and can provide 
brain activation maps with high spatial and temporal 
resolution. However, the basic limitations for each 

15 modality, such as the indirect nature of a f MRI measurement 
and the inverse problems for EEG and MEG, remain obstacles. 
In contrast, msMRI overcomes these limitations and directly 
measures magnetic sources originating from the neural 
firing with high spatio-temporal resolution. 

2 0 As discussed above, modeling has been performed that 

confirms the results described herein. A neuron in the 
brain typically consists of a single axon and multiple 
dendrites. Each dendrite or unmyelinated axon may be 
modeled using a current dipole model . In such a current 

2 5 dipole model, the current flowing inside a neuron is 
opposite to the current outside the neuron. The dipole 
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model considers effects of intracellular current ii and 
extracellular current i e and ignores the effect of through 
membrane current. It can be easily demonstrated that the 
through membrane current generates no net magnetic field 
5 outside neuronal membrane and can be ignored. For myelinated 
axons, the majority of neural current concentrates at nodes 
of Ranvier and is through membrane current. Contributions of 
a myelinated axon to overall magnetic field are small and 
can then be ignored. As the majority of axons in the brain 
10 are myelinated, dendrites are the main source of neuronal 
magnetic fields. 

Magnetic field generated by a current dipole at any 
observation point outside a neuron can be calculated by: 

15 where B is the strength of magnetic field, (i 0 is magnetic 
permeability in the space, r is the distance from the center 
of the dipole, p is current dipole, x denotes a cross 
product, and the bold font indicates a vector. 

The current dipole p points along the neuron in the 

2 0 direction of the advancing depolarization wave and has 
magnitude 

p = TttfcTiAVi ( 2 ) 
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where a is radius of the dendrite, o± is electrical 
conductivity of axoplasma, and Avi is the change of membrane 
potential . 

Because the length of a dendrite is much longer then 
5 the radius of the dendrite, the magnetic field generated by 
a current dipole at any observation point inside a dendrite 
can be conveniently calculated according to Ampere's 
circuital law, which states that a line integral of B around 
a closed path I is equal to the integral of current density 
10 j through any surface s enclosed by the path: 

jB*(B = t4> JJj*<fe (3) 

where • denotes a dot product. Considering the cylindrical 
symmetry of a dendrite and assuming a uniform current 
density, the magnitude of the field at observation point r 
15 is related to the intracellular current, it, by: 

LTia 

By combining Equations 1 and 4, the neuronal magnetic 
fields at any observation point can be calculated. 

In a simple two-dipole system with parallel (two 
20 dipoles in the same direction) and anti -parallel (two 
dipoles in the opposite direction) configurations, the 
magnetic field at any observation point is a vector 
summation of the magnetic field generated by each dipole at 
that point. Such a configuration is shown in FIG. 9. 
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Referring to FIG. 9, shown are configurations of a two- 
dipole system having parallel (left) and anti-parallel 
(right) configurations. Parameters used for calculating 
the magnetic fields are: radius of the dendrites a = 0.5 \im, 
electrical conductivity of axoplasma <j± = 2 Q^rrf 1 , and 

membrane potential Av± = 75 mV. 

As calculated, neuronal magnetic fields increase with 
distance r; peak on the surface of the dendrites; and then 
quickly decline. The fields are strongly localized, 
concentrated in and around the dendrites. 

Orientations of neuronal magnetic field, B n , may take 
any possible direction. The component of B n parallel to the 
B 0 (of the MR I scanner) will cause the spins at a point (x, 
y, z) in the transverse plane to acquire additional phases, 
(pi(x, y, z) , which depend on strengths of the local neuronal 
magnetic field, B n (x, y, z, t) . 

ft(x,j/,z) = j^B n// (x,y,z,t)dt= ^ yB nl (x, y, z, t) cos (0)dt (5) 

where TE represents echo- time, B n //(x, y, z, t) is the 
parallel component of B n (x, y, z, t) , 0 is angle between B 0 
and B n/ and y represents gyromagnetic ratio. 

The majority of Bni, which is the component of B n 
perpendicular to the B 0 , will have no net effects on the 
spins. Only a small fraction of Bni at the Larmor frequency 
will act like a Bi field to rotate the spins at a point (x, 
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y, z) away from the x-y plane. Additional phases, cp 2 (x, y, 
z) , generated by the neuronal magnetic field may be 
approximated by: 



5 where f is the fraction of Bni at the Larmor frequency. The 
value of f is typically very small. The <p 2 , therefore, can 
be ignored. 

Neuronal magnetic fields could be detected, at least 
in theory, by mapping either phases or magnitudes of MR I 

10 signals. The phase shift for a voxel is an integral of 
<pi(x, Yt z) over the voxel. It is interesting to note that 
the magnetic fields, B nr in the parallel dipole 
configuration are 180° out of phase. The average phase 
shift for an image voxel is then close to zero. Similarly, 

15 average phase shift for anti-parallel configuration is also 
near zero. 

An image voxel typically consists of many neurons. In 
a two adjacent two-dipole system in anti-parallel 
configuration, as shown in FIG. 10, magnetic fields in a 

20 two-dipole system are always opposite to that in the 
adjacent two-dipole systems. The integration of B n (average 
phase shift) over an image voxel is again approximately 
zero. While the analysis is based on a simple two-dipole 
model, the conclusion holds in more complex situations. The 

25 direction of magnetic fields on the left side of a current 




(6) 
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dipole is always opposite to that on the right. Thus, phases 
of MR I signals will always destructively add. In short, 
detecting neuronal magnetic fields by mapping MR I phase 
shifts will not succeed. 

Mapping neuronal activity with MR I by detecting 
magnetic fields induced by neural firing is 
straightforward. The magnitude of the MR I signal observed 
at a point (x, y, z) after 90° radio -frequency excitation is 
related to phase changes induced by the neuronal magnetic 
field by: 

s(x,y,z)cc p(x,y,z)cos(<p } (x,y,z)) ( 7 ) 

where p (x ,y , z) is spin density, and (pi(x, y, z) is defined 
by equations (5) . The MR I signal for an image voxel is an 
integral of s (x, y, z) over the voxel : 

S = £ £ y j1(x,y 9 z)dxdydz ( 8 ) 

where Ax, Ay, and Az is the dimensions of the voxel. By 
substituting Equations 5 and 7 to Equation 8, the MR I 
signal is related to TE by: 

Soz J J fp(x 9 y,z)cos( jjS tt (x 9 y 9 z 9 t)cos(0)dt)dxdydz (9) 

0 0 0 0 

Equation 9 implies that the MR I signal depends on both 
the strength and orientation of neuronal magnetic fields. 
When B n is perpendicular to B 0 , no MR I change will be 
induced by neuronal magnetic fields. 
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Because the component of B n perpendicular to B 0 does 
not generate significant MR I signal changes, focus may be 
placed on estimating the effects of the B n component 
parallel to the B 0 field. 
5 As discussed, an image voxel typically consists of 

millions of dendrites. Let us first consider that the 
dendrites are uniformly distributed in a voxel with 
parallel or ant i -parallel configurations and then 
generalizing the model later. 

10 MR I signal changes can be calculated by substituting 

the magnetic field distribution to Equation 9. When the 
dipoles are perpendicular to the B 0 field, the relationship 
between MR I signal changes and the number of dendrites 
firing simultaneously is plotted in FIG. 11. When the 

15 dipoles are parallel to the B 0 field, majorities of neuronal 
magnetic fields are perpendicular to the B 0 field. MR I 
signal changes are very small (near to 0) . 

The MR I signal change depends on orientations, 
configurations, and number of dendrites firing 

2 0 simultaneously. As used herein "orientation" means the 
spatial relationship between neuronal magnetic fields and 
the B 0 field. As used herein, "configuration" is the 
spatial relationship between dendrites (current dipoles) . 
When angles between dendrites change, it will change 

25 magnitudes and orientations of magnetic fields. As shown 
in FIG. 11, an anti-parallel configuration creates higher 
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signal changes compared to a parallel configuration. The 
number of dendrites firing simultaneously also plays an 
important role on the MR I signal changes. However, the 
exact number of dendrites firing at any given time during 
5 echo time TE is unknown and may be estimated. For FIG. 11, 
the parameters used for calculating the MR I signal changes 
are: radius of the dendrites a = 0.5 jam, electrical 
conductivity of axoplasma Oi =2 Q^m" 1 , TE = 100 ms, and 
membrane potential Avi = 75 mV, uniform spin density across 

10 the voxel. The number of dendrites here means the average 
number of dendrites fired during echo time TE for a typical 
MR I voxel. A conventional estimate is about 100,000 
pyramidal cells per square millimeter of cortex, each with 
tens to thousands of synapses. For a typical MR I voxel of 

15 3X3 mm 2 , synchronous activation of 0.1% to 1% dendrites 
could bring up the number of dendrites firing 
simultaneously to about 0.1 to 1 million. If so, the MR I 
signal changes should be about 0.44% to 3.6% for parallel 
configuration and 0.46% to 5.3% for anti-parallel 

20 configuration, which should be detectable. 

The above analyses are based on the assumption of 
rather simple configurations of dendrites: parallel or 
anti-parallel configurations. The real configurations of 
dendrites in the brain are typically much more complex. If 

25 the dipoles of a parallel configuration rotated by an angle 
a are still perpendicular to the B 0 field, MR I signal 
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changes are in the range defined by the parallel (the lower 
limit) and the ant i -parallel configurations (the upper 
limit) (FIG. 11) . However, if the dipoles being rotated 
are no longer perpendicular to the B 0 field, MR I signal 
5 changes will decrease. In the worst case where the rotated 
dipoles are parallel to the B 0 field, MR I signal changes 
will decrease by about 50%. 

So far, we have only considered a uniform distribution 
of dendrites in a voxel. The distributions of dendrites 

10 are actually not uniform in the brain. A general approach 
taken for the non-uniform distribution is sub-dividing a 
voxel into multiple subvoxels. Dendrite distribution in 
each subvoxel is approximated by a uniform distribution. 
The MR I signal change for the voxel can then be computed by 

15 calculating the weighted-average (weighted by the sizes of 
subvoxels) of the signal changes for each subvoxel. Our 
findings indicate that MR I signal change is dependent on 
the average dendrite density, but independent of dendrite 
distributions for the ant i -parallel configuration. For 

20 parallel configuration, the MRI signal change is 
insensitive to dendrite distributions for low dendrite 
densities (<15000 dendrites/mm 3 ) and are slightly smaller 
than that for uniform distribution when dendrite densities 
are high (>15000 dendrites/mm 3 ) . 

25 Magnetic fields generated by neural firing on the 

scalp detectable by MEG based on the model may be 
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estimated. The magnetic fields outside the brain are 
strongly dependent on the configurations of dendrites. An 
anti-parallel configuration will generate no net magnetic 
field on the scalp. A parallel configuration will generate 
a magnetic field of 7 x 10" 12 T on the scalp 4 cm away from 
the current source, assuming no attenuation by the scalp 
and one million dendrites firing simultaneously. The 
orientations of dendrites in the brain are very complex and 
are neither simply parallel nor simply anti-parallel. The 
magnetic fields on the scalp should then be somewhere 
between that generated by an anti-parallel configuration 
and that generated by a parallel configuration, i.e., 
between 0 to 7 x 10" 12 T. MEG detection is also dependent 
on the orientation of neuronal magnetic fields relative to 
the MEG detector. When the neuronal magnetic field is 
parallel to the detecting coil, no signal will be detected. 
The actual magnetic field reported by a MEG experiment is 
about 10" 13 T and is in the range estimated based on our 
model . 

For BOLD contrast or other magnetic field 
inhomogeneity, the effects on MR I signals are typically 
modeled by a T 2 * effect. The relative signal change created 
by a BOLD contrast is described by: 

AS = ^y^*R{TE (10) 
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where S 0 and S x are MR I signals without and with BOLD 
effect, and R' 2 is the external decay rate due to the BOLD 
contribution. The R' 2 is independent to TE . 

The TE dependence of MR I signal changes due to neuronal 
5 magnetic fields can be computed based on Equation 9. The 
relationship of TE and MR I signal changes is given by 
Equation 11 and is plotted in FIG. 12, which shows computed 
MR I signal changes with TE for parllel (P) and ant i -parallel 
(AP) configurations. It is very interesting to see a 
10 nonlinear relationship between TE and MR I signal changes. 
When TE doubled, the MR I signal changes increase by an 
average 3.8 times. The unusually strong dependence of MR I 
signal changes on the echo-time TE provides an excellent 
opportunity to experimentally validate the model. 

AjcAvAz TE 

\ I z) cos[ \yB n (x ,y,z, t) cos{0)dt]dxdydz 

15 AS = ^L = ±JU> _J 1 (11) 

J J jp(x 9 y,z)dxdydz 

0 0 0 

To verify the TE depedence, experiments were performed 
on three subjects for three different TEs . The signal 
changes for all activated brain areas were measured. The 
experimental data show a nonlinear relationship between TE 
2 0 and MR I signal changes (FIG. 8) and are supportive to the 
theoretical modelling . 

These results demonstrate that MR I signal changes 
induced by neuronal magnetic fields are in detectable 
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levels. A MR I magnitude map may provide a better result 
than a phase map. Amplitudes of msMRI signals depend on 
orientations, configurations, and density of dendrites. 
Furthermore, the MR I signal changes have a nonlinear 
5 relationship with TE. The strong dependence of MR I signal 
changes on TE have been confirmed by experimental results. 

While the present invention has been described with 
respect to a limited number of embodiments, those skilled 
in the art will appreciate numerous modifications and 
10 variations therefrom. It is intended that the appended 

claims cover all such modifications and variations as fall 
within the true spirit and scope of this present invention. 
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